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We study theoretically the nonlinear four-wave mixing response of an ensemble of coupled pairs of 
quantum dots (quantum dot molecules). We discuss the shape of the echo signal depending on the 
parameters of the ensemble: the statistics of transition energies and the degree of size correlations 
between the dots forming the molecules. 

I. INTRODUCTION 

Optical experiments on quantum dot molecules (QDMs) 0, H, Q have shown that QDMs differ in many ways from 
single quantum dots (QDs). In a recent work [1| we have discussed the effects of collective interaction between closely 
spaced dots and the modes of electromagnetic field (sub- and superradiance @), as manifested in the luminescence 
decay from a single QDM. Here we will consider the nonlinear four-wave mixing (FWM) response from an ensemble 
of QDMs 0. 

The FWM spectroscopy is a powerful tool for extracting information on lifetimes and homogeneous dephasing from 
inhomogeneous ensembles [f| Q. In the case of Markovian decoherence, the interpretation of the measured FWM 
signal is straightforward [1] . More care must be taken for non-Markovian pure dephasing, where the FWM response 
does not reproduce the optical response of a single system [9| . As we shall see, the time evolution of the FWM signal 
from QDMs is much more complicated than that from QDs and carries information on the interaction between the 
dots and on the correlations of their size. Therefore, before one can proceed to the discussion of decoherence effects, 
it is necessary to analyze all the factors that affect the shape of the FWM response. This is the goal of the present 
paper. We will study the shape of the FWM signal for a pair of QDs and its dependence on the size distribution and 
size correlations between the dots, as well as on the interaction between them, in absence of any decoherence. 

II. THE MODEL 

Each QDM will be modelled as a four-level system with the basis states |00), |01), |10), |11), corresponding to the 
ground state (empty dots), an exciton in the second and first QD, and excitons in both QDs, respectively. A FWM 
experiment consists in exciting an ensemble of QDMs with two ultrashort laser pulses, arriving at t\ = — r and ti = 0. 
We assume that the pulses are spectrally very broad to assure resonance with all the QDs in the ensemble. First, 
consider a single QDM composed of two QDs with energies = E ± A. We will describe its evolution in a frame 
rotating with the frequency E/H. Then, in the rotating wave approximation, the Hamiltonian for this single QDM is 

H = A(|l)(l|®I-I®|l>(l|) + y(|01)(|10| + |10)(|01|) (1) 

+ \ E [e-^' +EU) (|0)(1| ® 1 + 1 ® |0)<1|) + H.c. 

i 

where fi and §i are the amplitude envelopes and phases of the pulses, respectively, V is the coupling between the 
dots (tunnel or Forster), I is the unit operator and the tensor product decomposition refers to the two QDs making 
up the molecule. The pulses do not overlap in time. 

III. THE SYSTEM EVOLUTION AND FWM RESPONSE 

If the durations of the pulses are much shorter than both ft/A and h/V, the action of each of them corresponds 
to the independent rotation of the state of each QD, that is, to the unitary transformation U; = Ui <g> Ui, where 
Ui = cos(aj/2)E — ism(ai/2)[e~ t ^ i+Et ^\0){l\ + H.c] and a t — dtfiit) is the pulse area. Between the pulses and 
after the second pulse the system evolution is generated by the time-independent Hamiltonian given by the first two 
terms in Eq. |T]). The resulting evolution operator reads 

W(i) = |00)(00| + |ll)(ll|+cos(Ai)(|01)(01| + |10)(10|) 

+isin(Ai) [cos(20)(|Ol)(Ol| - |10)(10|) - sin(20)(|Ol)(lO| + |10)(01|)], 



where sin(20) = V/VV 2 + A 2 . 
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The system state at a time t > is represented by the density matrix 

p(t) = W(t)U 2 W(r)Ui jO (0)[W(t)U 2 W(T)Ui] t , 

where p(0) — |00)(00| is the initial state. The optical polarization of a single QDM under consideration is proportional 
to P{t) = (ll|p(t)(|01) + 1 10)) + (<01 1 + (10})/?(t)|00) + c.c. In order to extract the FWM polarization we pick out 
the terms containing the phase factor e^ 2 ^ 2 ^ 1 ) . The total optical response from the sample is obtained by summing 
up the contributions from individual QDMs with a weight factor g(Ex,E2) reflecting the distribution of transition 
energies in the sample (we neglect a possible variation of dipole moments). As a result we obtain the third-order 
FWM signal P FWM (t) = ^ P>(t) + c.c, where 



Pi(t) 

p 4 (t) 



isinaisin 2 y J dE 1 dE 2 g(E 1 , E 2 )e~ lE( - f ~ T ) cos [Q(t - r)] , 

— sin 2ai sin 2 a 2 J dE 1 dE 2 g(E 1 , E 2 ) e - lE{t - r) sin [Q(t - r)] , 

^sin2aisin 2 y cosa 2 / dE\dE 2 g{E\, E 2 )e~ iE ( t ~ T " 1 sin [fi(t + r)], 

i^-sinaisin 2 ^cosa 2 € tdE 1 dE 2 g(E 1 , E 2 ) e - zE{t - T) 
il z 2 

cos [fi(f - r)] + cos [fi(f + t)] 



where = yV^ 



A 2 . We assume a Gaussian distribution function 

g{E 1 ,E 2 ) = 
1 

: exp 



(2a) 
(2b) 
(2c) 

(2d) 
(3) 



27R7 2 V / T 3 P 2 



Pi - Pi) 2 - 2p (Pi - Pa) (P 2 - P 2 ) + (P 2 - P 2 ) 



2(l-p 2 )a 2 



with identical energy variances for both QDs (consistent with the symmetric photoluminescence spectrum observed 
in the experiment [l|) and a correlation coefficient p. Note that this distribution corresponds to an uncorrelated 
Gaussian distribution of the parameters P and A with variances <je = a°J (1 + /°)/v / 2 and cta = <J\J (1 — p)/y2 (a 
small variance of the energy difference A implies correlated energies). For the calculations we will use the values for 
the energy variance a — 8 meV and for the average energy mismatch A = (Pj — P 2 )/2 = 11 meV [![. In order to 
fix the relative amplitudes of the contributions Pi{t) we assume a, <C 1 and expand the trigonometric coefficients in 
Eqs. ((Md). 

The detection of weak signals originating from QDs is based on a heterodyne technique [6| : The response Pfwm is 
superposed onto a reference pulse E Te f(t) = f re f(t — to) e ~ 1 ' E( ' t ~* ' > + c.c, arriving at a time to. We assume a Gaussian 
envelope f T ei(t) — exp[— (l/2)(t/r rc f) 2 ]/(v / 27rT ro f). The measured signal is proportional to 



P(to,r) 



dtP£> M (t)E£(t) 



(4) 



where PpwM an d ^ref^ are ^ ne positive frequency part of the FWM signal and the negative frequency part of the 
reference pulse, respectively. 

In the case of decoupled QDs (V = 0), the only non- vanishing contribution to Pfwm is Pi- Inserting Eq. © into 
Eq. (|2a|) and using Eq. (|4]) one finds 



n(-) 



F{t ,T) 



: exp 



q 2 (T-* ) 2 + T 2 ef A 2 
2(^> 2 + l) 



A(r - 1 ) 



'rof L 



This result is plotted in Fig.[T^,,d for an ultrashort reference pulse and for a pulse with a realistic duration of r rc f = 43 fs, 
i.e., 100 fs full width at half maximum (FWHM). The obtained result is remarkable for two reasons. First, the shape 
of the echo peak for noninteracting dots does not depend on the correlations between the transition energies Pi and 
P 2 . Second, the beats related to the energy mismatch are always in phase with the FWM echo, so that they do not 
lead to oscillations in the time-integrated signal. The only trace of these beats are the oscillations in the tails of the 
echo peak. Third, the amplitude of the echo contains the factor exp[— (l/2)r r 2 f A 2 /(t 2 f er 2 + 1)] and becomes small if 
the QDM is formed by very different dots. 
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FIG. 1: FWM signal from the QDM sample detected with an ultrashort reference pulse (a-c) or with a 100 fs FWHM pulse 
(d-f). Solid lines: p = 0; dashed lines: p — 0.875. 
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FIG. 2: Time-integrated FWM signal as a function of the delay for two values of the coupling as shown and for p = (solid) 
and p = 0.875 (dashed). 



If the QDs are coupled, the other three terms contribute to the signal. This results in a modification of the shape 
of the echo peak, which now depends on the degree of correlation between the energy parameters of the two dots. In 
Fig. [T] we show the results obtained by numerical integration. For V = 0.8 meV (Fig. []J>,e), which is a realistic value 
of the Forster coupling for 5 nm distance between the dots, the shape is only slightly modified. A much larger effect 
is observed for stronger coupling (as might result from tunneling). For the parameters chosen here, the characteristic 
sharp feature appearing at the top of the response peak for strong couplings is considerably smeared out by the 100 
fs reference pulse (Fig. [TJ;,f). Nevertheless, the evolution of the echo shape with growing coupling remains clear. 

As mentioned above, in the limit of vanishing interaction there are no beats in the time-integrated FWM signal as 
a function of the delay r. This changes for V ^ due to the contribution from the term P-$(t) and second part of 
P±(t) (Fig. [2). At the center of the echo peak, t = r, these terms have the variable phase factors cos2Slr and sin2f2r. 
As a result, they lead to oscillations of the integrated FWM signal as a function of r but only as long as tcta ^ 1- 
For longer times these terms do not contribute because of the phase averaging due to inhomogeneous distribution of 
energy differences A. Such oscillations of the time-integrated signal might be a useful signature of interaction between 
the QDs. Unfortunately, for the system parameters used here they are present only for r ^ 1 ps and will be covered 
by phonon-related features |lj. 



IV. CONCLUSION 



Our results show that the shape of the time-resolved FWM signal may provide some clue on the properties of the 
QDMs in the ensemble, including coupling and, for coupled dots, size correlations between the QDs. We have studied 
only the amplitude of the FWM response. It should be noted, however, that the terms P2 and P3 are real, while the 
other two are purely imaginary. Therefore, much more information is contained in the phase properties of the signal. 
Moreover, nonlinear experiments can be performed beyond the perturbative regime [10] . This might allow one to 
selectively investigate selected contributions to the FWM signal by exploiting their different dependence on the pulse 
areas oti^- 
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